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Abstract Turbulence in an incompressible fluid with and without a magnetic field as well as moder- 
ately compressible MHD turbulence are compared. The results of three numerical simulation models 
in three dimensions of resolution up to 512 3 are used for this purpose. The compatibility of the spec- 
tra of all three models with the Kolmogorov spectrum is confirmed. For the magnctohydrodynamic 
(MHD) models the probability distribution functions of the velocity components perpendicular to the 
external magnetic field are like the incompressible hydrodynamic (HD) model while those parallel to 
the field have a smaller range of velocities. The probability distribution functions of the transverse 
velocity increments for the MHD models decline slower than the incompressible HD model. The simi- 
larity of incompressible HD and both incompressible and compressible MHD turbulence persists over 
high order longitudinal structure function scaling exponents measured in the global reference frame as 
well as for motions perpendicular to the local mean field. In these two frames the longitudinal scaling 
exponents of both MHD models seem to follow theoretical incompressible HD dissipation structure 
predictions while the transverse scaling exponents of the incompressible MHD model seem to follow 
the predictions for incompressible MHD. In the local magnetic system the motions parallel to the local 
mean field for both MHD models are different from incompressible HD motions. The fields of the MHD 
simulations are decomposed into Alfvenic, fast, and slow modes. In the global reference frame and for 
motions perpendicular to the local mean field the Alfvenic mode is mostly responsible for the fact 
that the longitudinal components of both MHD models follow the incompressible HD model and the 
incompressible HD theoretical dissipation structure predictions. In the global reference frame for the 
incompressible MHD model both the Alfvenic and slow modes seem to contribute significantly to the 
fact that the transverse components follow the incompressible MHD theoretical predictions. For mo- 
tions perpendicular to the magnetic field for the incompressible MHD model it appears that the slow 
mode contributes most to the fact that the transverse component seems to follow the incompressible 
MHD predictions. 

PACS 47.27.E- • 52.30.Cv • 52.35.B • 95.30.Qd 



1 Introduction 

Given that turbulent motions (velocities) are inherently stochastic, velocity field statistics are an essen- 
tial tool for extracting its physical properties. In 1941 Kolmogorov [lj published his theory (hereafter 
K41) on the local velocity field structure of turbulence for incompressible hydrodynamic (hereafter 
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IHD) fluids with high Reynolds numbers, in which he hypothesized the existence of homogeneity, 
isotropy, and an inertial range (hereafter K41). The inertial range in the kinetic energy spectrum is 
a range in wave number, k ~ l/l where I is the separation length between two points often regarded 
as the eddie size, which lies between the typical large scales and the dissipative scales and is there- 
fore independent of external input such as turbulence drive or output such as dissipation [2]. In this 
range the spectrum is also independent of viscosity and exhibits a power-law dependence of the form 
E(k)=Ke~ 2 / 3 /c~ 5 / 3 , where K is the universal Kolmogorov constant and e is the average rate of energy 
dissipation per unit mass. 

It is known that the inclusion of an external magnetic field in the turbulence of an incompressible 
fluid brings about many substantial complications and can drastically modify the turbulent motions. 
The inclusion of compressibility also causes significant changes in the properties of turbulence. Gotoh, 
Fukayama, & Nakano [3] thoroughly and successfully investigated the velocity field statistics of high- 
resolution direct numerical simulations for IHD turbulence. The purpose of this paper is to address 
the question, are the velocity field statistics completely different from IHD when a magnetic field and 
compressibility arc included? We extend the analysis of Gotoh et al. [3] to incompressible magnetohy- 
drodynamic (hereafter IMHD) and compressible magnetohydrodynamic (hereafter CMHD) turbulence 
in an attempt to answer this question. This is considered mostly through the use of high order scaling 
exponents of the velocity structure functions. 

Emphasis is also placed on the differences between the scaling exponents of longitudinal and trans- 
verse structure functions. For high order scaling exponents the findings of Gotoh et al. [3] are confirmed 
for structure functions measured in the global reference frame for IHD turbulence. Scaling exponents 
are then explored for IMHD and CMHD turbulence in the global frame and the local frame, defined by 
the local mean magnetic field. The scaling exponents are then decomposed into Alfvenic, fast, and slow 
modes and the contributions of individual modes are explored. The paper is organized as follows. Sec. 2 
presents the numerical aspects of the three numerical simulation models considered (IHD, IMHD, and 
CMHD). Sec. 3 contains analysis of the kinetic energy spectra of the three models. In Sec. 4 probability 
distribution functions are considered. Structure functions are explored in Sec. 5. Sec. 6 presents scaling 
exponents of the structure functions. In Sec. 7 the scaling exponents are decomposed into the three 
MHD modes. Finally, Sec. 8 includes a discussion and summary. 



2 Simulation Specifics 

The IHD data used in this study came from simulations with a grid size of 256 3 . For both MHD 
models, a high resolution grid size of 512 3 was used and the external magnetic field is in the positive 
x-direction. 



2.1 Incompressible HD and MHD 

The IMHD code described below is used for the IHD model with B set to zero. We have calculated the 
time evolution of incompressible magnetic turbulence subject to a random driving force per unit mass. 
We have adopted a pseudospectral code to solve the incompressible MHD equations in a periodic box 
of size 2ir: 

n 

— = (V x v) x v - (V x B) x B + t/V 2 v + f + VP', (1) 

^=Vx (vxB)+7?V 2 B, (2) 

V ■ v = V ■ B = 0, (3) 

where f is a random driving force, P' = P/p + v • v/2, v is the velocity, and B is magnetic field 
divided by (Airp) 1 / 2 . In this representation, v can be viewed as the velocity measured in units of the 
r.m.s. velocity, v, of the system and B as the Alfven speed in the same units. The time t is in units 
of the large eddy turnover time (~ L/v) and the length in units of L, the inverse wavenumber of 
the fundamental box mode. In this system of units, the viscosity v and magnetic diffusivity r\ are the 
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inverse of the kinetic and magnetic Reynolds numbers respectively. The magnetic field consists of the 
uniform background field and a fluctuating field: B = Bo + b. We use 21 forcing components with 
2 < k < V^12, where wavenumber k is in units of L^ 1 . Each forcing component has correlation time 
of one. The peak of energy injection occurs at k « 2.5. The amplitudes of the forcing components are 
tuned to ensure v ~ 1. We use exactly the same forcing terms for all simulations. The Alfven velocity 
of the uniform background field, B , is set to 1. We consider only cases where viscosity is equal to 
magnetic diffusivity: 

v = ij. (4) 

In pseudo spectral methods, the temporal evolution of equations |T]) and ^ are followed in Fourier 
space. To obtain the Fourier components of nonlinear terms, we first calculate them in real space, 
and transform back into Fourier space. The average kinetic helicity in these simulations is not zero. 
However, previous tests have shown that our results are insensitive to the value of the kinetic helicity. 
In incompressible fluid, P' is not an independent variable. We use an appropriate projection operator 
to calculate VP' term in Fourier space and also to enforce divergence-free condition (V-v — V-B = 0). 
We use up to 256 3 collocation points. We use an integration factor technique for kinetic and magnetic 
dissipation terms and a leap-frog method for nonlinear terms. We eliminate the 2At oscillation of the 
leap-frog method by using an appropriate average. At t = 0, the magnetic field has only its uniform 
component and the velocity field is restricted to the range 2 < k < 4 in wavevector space. 

Hyperviscosity and hyperdiffusivity are used for the dissipation terms. The power of hyperviscosity 
is set to 6, so that the dissipation term in the above equation is replaced with 

-^ 6 (V 2 ) 6 v, (5) 

where i/g is determined from the condition Vh(N/2) 2h At ss 0.5 (see [4]). Here At is the time step and iV 
is the number of grid points in each direction. The same expression is used for the magnetic dissipation 
term. Cho & Lazarian [5] can be referenced for more details on the above discussion. 



2.2 Compressible MHD 

We use a third-order accurate hybrid essentially non-oscillatory (ENO) scheme (see [6]) to solve the 
ideal isothermal MHD equations in a periodic box: 

dp/dt + V • (pv) = 0, (6) 
dv/dt + v Vv + p-^i^p) - (V x B) x B/47rp = f , (7) 

dB/dt - V x (v x B) =0, (8) 

with V • B = and an isothermal equation of state. Here f is a random large-scale driving force, 
p is density, v is the velocity, and B is magnetic field. The rms velocity SV is maintained to be 
approximately unity (in fact SV ~ 0.7), so that v can be viewed as the velocity measured in units 
of the r.m.s. velocity of the system and B/^/47rp as the Alfven velocity in the same units. The time 
t is in units of the large eddy turnover time (~ L/SV) and the length in units of L, the scale of the 
energy injection. The magnetic field consists of the uniform background field and a fluctuating field: 
B = Bo -l-b. For our calculations we assume that Bo/^Airp ~ SB/y/Anp ~ SV. The sound speed is one 
therefore the resulting Mach number and plasma (5 are of order unity. As is the case with other finite- 
difference numerical schemes, our compressible MHD code inherits numerical viscosity (and diffusion), 
which results in energy dissipation on small scales. See [7] for details on the above discussion. 



3 Kinetic Energy Spectrum 

The three-dimensional kinetic energy spectrum for each model, taken after turbulence has reached 
a statistically stationary state, is shown in Figure [T] The inertial range of the spectrum is of most 
interest for this study. The shape of the curves at wave numbers less than those of the inertial range 
depends on the energy at the injection scale and how the turbulent motions are excited. The shape 
of the curves at wave numbers greater than those of the inertial range is affected significantly by 
viscosity, the bottleneck effect [5] , and the resolution of the model. Despite the fact that the Kolmogorov 
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phenomenology, E(k)cxfc~ 5 / 3 , is derived for IHD, in the inertial range all three curves demonstrate a 
very similar power law dependence, suggesting that the energy transfer process is consistent over this 
range of wave numbers. This may be surprising considering the many changes resulting from the 
inclusion of a magnetic field, such as an Alfvenic mode propagating parallel to the magnetic field and 
elongation of the structures in the velocity field. The inclusion of compressibility introduces significant 
changes as well, such as shocks. Goldreich & Sridhar [9] provide an explanation of this behavior. See 
[TO] also for a review. Deviations from the 5/3 law have become an issue of debate (see [TTUT2"] . and 
references therein). We do not address this issue here. 



4 Probability Distribution Functions 

The probability distribution function (hereafter PDF) is used to obtain information about one-point 
velocity statistics (see |19j for detailed discussions on astrophysical applications). The single-point 
PDF's of different components of the velocity field for each model shown in Figure [2] are histograms, 
normalized to the grid size, of the values of the velocity at each point. For the IHD model the x, 
y, and z components were averaged. For the two MHD models, the histograms were separated into 
components parallel (x direction) and perpendicular (y & z directions) to the external magnetic field, 
(B ex t); therefore, the PDF's for the x components are presented separately and the y and z components 
are averaged. 

The PDF for the IHD model is somewhat asymmetric. Any asymmetry or skewness in the PDF's 
should disappear when averaged over multiple time steps [3J. The PDF's of the MHD simulations 
parallel to B ext are similar to each other, fairly symmetric, wider toward the top than a gaussian, and 
at large amplitudes decay faster than a gaussian. Their slopes are the steepest of the PDF's presented 
and given that the area under the curve must be the same, their ranges of velocities are smaller. For 
IMHD turbulence the divergence of the velocity vector is zero meaning that there are no sources, sinks, 
or shocks and that only eddies can exist. This is not the case for CMHD simulations in which shocks 
are produced at high velocities which dissipate the energy perhaps contributing to the PDF of the 
parallel CMHD model having the least high velocity points. 

The PDF's of the IMHD and CMHD models perpendicular to B ext approximately follow each other 
and are even more similar to each other than the parallel components. They are very symmetric and 
gaussian-like, have the most shallow slopes of all, and therefore have the most high velocity points. 

When an external magnetic field is introduced into the turbulent motions, an Alfvenic mode is 
produced with perturbations perpendicular to the direction of the field. This additional contribution 
to the velocity is responsible for the higher velocities found in the perpendicular PDF's. The additional 
kinetic energy given to the perpendicular direction explains the separation of these two components 
seen in the top of Figure O This appears to be equally valid for both positive and negative values of 
the velocity. 

An average over 10 time steps of the velocity PDF's for a CMHD simulation with 256 resolution 
was performed. The 10 parallel time steps were averaged as well as the 20 perpendicular time steps, 
10 for y and 10 for z. These PDF's are plotted in Figure [3J along with the CMHD PDF's mentioned 
above. The PDF's are normalized to match their maxima in order to better compare their tails. The 
parallel and perpendicular PDF's for the single time step lie sufficiently close to the expected random 
errors obtained from the PDF's averaged over multiple time steps to justify using a single time step 
and averaging the perpendicular components. 

Figures [4] and [5] show PDF's of the transverse increments of the velocity field which are histograms, 
normalized to the grid size, of the values of the difference between the velocity field at each point and 
the field at an adjacent point along the direction specified. Positive values suggests that v is increasing 
in the positive direction while negative values suggest that it is decreasing. In Figure [4] the increments 
of the x-component of the velocity field taken along the y and z directions, parallel to B ext for the MHD 
models, are averaged. Figure O shows the increments of the y-component along the x and z directions. 
The results are not averaged given that one is parallel to B ext while the other is perpendicular. 

In both Figures, all three models demonstrate near perfect symmetry. This is true despite the 
orientation with respect to B ext . In Figure [5] as expected for the IHD model, x and z are nearly 
identical because there is no magnetic field to bring about differences. On the other hand for the MHD 
models, x extends to higher values of the velocity increment than z. 
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Because of the fact that the difference in range between the x components for the IMHD and 
CMHD models is much greater than the difference in range between the z components for the IMHD 
and CMHD models, it is possible that most of the increase in range in going from IMHD models to 
CMHD models by introducing compressibility is coming from the components parallel to B ext . This is 
supported by the fact that for the PDF's in Figure [2] the perpendicular components follow each other 
closely while the parallel components are quite different from each other. 

5 Structure Functions 

Structure functions are used to obtain information about two-point velocity statistics. Structure func- 
tions describe the distribution of velocity structures of size I and are defined as 



where Sv is a velocity increment and (...) represents ensemble averaging. Structure functions are defined 
with the absolute value of the velocity increment to ensure their positivity for all orders. Structure 
functions are separated into a longitudinal component naturally defined by the velocity increment 



and transverse components Svt{1) chosen randomly in the plane perpendicular to the longitudinal 
component. When higher order structure functions are used the PDF tails contribute more, where as 
you can see in Figure [3] the errors are much larger. 

For each separation length I from 1 to 128, Sv is calculated at 100,000 points randomly distributed 
throughout the entire periodic velocity field. In the global reference frame the direction of I, or the 
direction of the longitudinal structure function, is also random. The structure functions for the IHD 
model are calculated in the global reference frame. For the MHD models I can be chosen randomly 
(global frame), parallel to the local mean magnetic field, or perpendicular to the local mean field. 

It is well known that, for MHD turbulence in a strongly magnetized medium energy cascade occurs 
mainly in the direction perpendicular to the mean field in the Fourier space, which makes eddies 
elongated [T3 l fr4 l H5 l [T6] along the global mean magnetic field. However, according to [17] and [6], for 
MHD turbulence in a strongly magnetized medium eddies are not aligned along the direction of the 
total mean field but rather are aligned along the local mean field. This being the case, the local mean 
field is the physically relevant background for the eddie dynamics. The local frame is separated into 
two types which are represented in Figure El In (a) the points from which the structure functions are 
calculated are chosen along the local mean field, this configuration will be referred to as local parallel. 
In (b) the points are chosen in the plane perpendicular to the local mean field, this configuration will 
be referred to as local perpendicular. 

Figure[7]shows the second order longitudinal, transverse, and total structure functions of the velocity 
field. The column on the left is the global frame, the middle column is the local parallel frame, and 
the column on the right is the local perpendicular frame. In the local frames the structure functions 
of the IHD model calculated in the global frame are shown for comparison. 

The inertial range of the turbulent motions manifests itself in the structure functions as the range 
of / over which the slope is constant in a log- log plot. The first maximum of a structure function is 
related to the size of the largest structures in the corresponding field, the radius of the largest eddies 
for example. It will be found at the separation length at which the largest Sv is observed. 

We confirm the result of [3] that in the global frame for IHD turbulence the values of the transverse 
structure functions are greater than those of the longitudinal. In extending our analysis to IMHD and 
CMHD we find that the structure functions of both MHD models in the global and local perpendicular 
reference frames exhibit this behavior as well. For each MHD model in the local parallel reference 
frame the values of the longitudinal and transverse structure functions are very similar. In the global 
and local parallel frames the IMHD structure functions are greater than those of the CMHD model 
at lower values of the separation length where they appear more like the IHD structure functions in 
the global frame. In summary, the MHD models in the global frame and in the plane perpendicular to 
the local mean field are significantly more similar to the IHD structure functions than the structure 
functions of the MHD models in the local parallel frame. Also, the shapes of the structure functions 
in the global and local parallel frames for the IMHD model are more like that of the IHD model than 
the CMHD model is. 



srW(i) = (\5vQ)\*>), 



(9) 



foi(0 = [v(x + I)-v(x)H, 



(10) 



6 



Hall et al. 



6 Scaling Exponents 

Self-similarity or scale-invariance implies that fluid turbulence at larger scales can be reproduced by 
the magnification of turbulent motions at smaller scales. At the dissipation scales self similarity fails 
as turbulence forms non-Gaussian dissipation structures [2J. Self-similarity is even an approximation 
in the inertial range where turbulent motions become increasingly sparse in space and time at smaller 
scales, a property known as intermittency. In 1962 Kolmogorov published a modification of his 1941 
theory including the effects of intermittency or scale- variance |18) . One way to investigate intermittency 
is by studying the scaling exponents of velocity fluctuations in the inertial range, given that for a scale- 
invariant flow the scaling exponents are linear functions of the order of the structure functions. 

Despite the fact that the origin is not yet fully understood, the scaling behavior of structure func- 
tions is very useful and interesting (see [19j for a comprehensive discussion on astrophysical situations; 
see also [2D])- Within the inertial range, a structure function of arbitrary order scales with the sepa- 
ration length I as S^(l) cx l^ p \ From this scaling relation and a few transformations self-similarity 
between two structure functions of arbitrary orders, p and p', can be obtained: 

SW(l) oc [S , (p')(0] c(p)/c(p ' ) , (11) 

demonstrating how the structure functions S^(l) and S^ p ^(l) are related (see [2j). It also brings to 
light the possibility of using a structure function of an arbitrary order p' as a normalization. In Figure 
[5] the sixth order longitudinal structure function is plotted as a function of the third order longitudinal 
structure function where I is parallel to the local mean magnetic field, demonstrating how self similarity 
can be used to obtain scaling exponents. A minimum chi-square linear fit was used to determine C(6). 
The scaling relation S^ p \l) oc fiw in the inertial range has been derived analytically in Kolmogorov 
phenomenology yielding ((p) = p/3. For this reason the third moment is used as the normalization. 
If the effects of dissipation on the structure functions are independent of the order one may conclude 
that not only the third order moment but any moment S^ p \l) oc ' can be used [2J. 

She & Leveque [22] proposed a scaling relation that contains three parameters (see also [2Tl[23lf24] ) : 
g related to the scaling Sv ~ l~5 , x related to the energy cascade rate t~ l ~ l~ x , and C the codimension 
of the dissipative structures: 

C(p) = ?(l- a! )+C[l-(l-^)f]. (12) 
For IHD turbulence using g = 3, x — |, and C — 2 She & Leveque [22] obtained 

C SL \P)=1 +2[1 -(§)«], (13) 

suggesting that dissipation occurs over one-dimensional structures, such as filaments or vortices (here- 
after SL1). For three-dimensional IMHD turbulence Muller & Biskamp [23] (see also [25] ) in extending 
the She & Leveque [22] model to IMHD turbulence and using g = 3, x = |, and (7 = 1 obtained 

C SL \p) = l + l-(^) f (14) 

assuming the dissipative structures to be two-dimensional, such as sheets (hereafter SL2). Finally, 
for CMHD turbulence Padoan et al. [26] showed that the velocity field for subsonic MHD turbulence 
follows SL1, while for highly supersonic MHD turbulence it follows SL2. 

Figures l9l through fTTI contain the scaling exponents as a function of the order of the longitudinal, 
transverse, and total structure functions normalized to the third order. For the MHD models the 
Figures are organized as follows. In Figure [9] the global reference frame is used, in Figure [TOl the local 
parallel reference frame is used, and in Figure [TT] the local perpendicular reference frame is used. For 
the IHD model the global scaling exponents are also presented in Figures [10] and [Tl] for comparison, 
given that the local reference frame has no meaning without a magnetic field. Despite the fact that 
the error in the scaling exponents increases with order from the increasing fluctuation in the structure 
functions, emphasis is placed on high order structure functions (p > 4) given that at low orders the 
theoretical predictions considered here differentiate themselves very little. 

Beginning with Figure [9] in the global reference frame the longitudinal scaling exponents of all 
three models follow SL1 quite well. In the transverse case although the slope of the IHD model initially 
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follows SL1 it begins to flatten at around p — 7. Its shape is in agreement with the results of [3]. 
The IMHD model follows SL2 quite well while the CMHD model is between the two. The IHD results 
mentioned here agree with those obtained by |3] in that for p > 4 the longitudinal scaling exponents 
fall slightly below SL1 while the transverse exponents are significantly lower. For all three models, 
despite the fact that the longitudinal scaling exponents follow SL1 much better than the transverse 
follow SL2, the scaling exponents calculated from the total structure functions are much more SL2-like 
than SLl-like for all three models. Their behavior is very similar to the transverse case, suggesting that 
the transverse scaling exponents have a significantly larger amplitude than the longitudinal, which can 
be seen to be true in Figure [7j The same holds for Figures [10] and [TT] 

In Figure [10] in the local parallel reference frame for the longitudinal component, parallel to the 
local mean magnetic field in this case, the IMHD model is between SL1 and SL2 while the CMHD 
model mostly follows SL2. For the transverse component, perpendicular to the local mean magnetic 
field in this case, the IMHD model is between SL1 and SL2 while the CMHD model follows SL2. 

In Figure [TT1 in the local perpendicular reference frame for the longitudinal component, perpen- 
dicular to the local mean magnetic field in this case, the best agreement to SL1 is found for all three 
models at all values of p. For the transverse component, also perpendicular to the local mean magnetic 
field in this case, we do not find a similar agreement; the IMHD model is more SL2-like while the 
CMHD model remains in between. 

We confirm the results of [3] that for IHD turbulence the transverse scaling exponents, (t(p) are 
less than the longitudinal Cl(p) for orders greater than four. In extending the analysis to incompress- 
ible and compressible MHD turbulence we find that in the global reference frame and in the local 
perpendicular reference frame this behavior continues. This is not necessarily the case in the local 
parallel reference frame. This difference is probably related to the fact that in the local parallel frame 
the longitudinal component necessarily probes the direction parallel to the local mean field while the 
transverse component probes the plane perpendicular to the local mean field, thus incorporating a 
consistent difference between the two components that does not exist in the other two frames. 

The most conspicuous SL1 cases above are both longitudinal, namely Cl°(p) an d Cl° C ^(p) m 
which all models follow SL1 suggesting physically that the dissipation structures are one dimensional 
hydrodynamic-like vortices. The most conspicuous SL2 cases are both the transverse cases correspond- 
ing to the longitudinal SL1 cases above, namely Ct°(p) ano - Ct C ± (p) m wm ch the IMHD model mostly 
follows SL2 suggesting physically that the dissipation structures are two dimensional sheets dominated 
by MHD processes. It is possible and perhaps more likely that the dissipation structures only appear 
one or two dimensional in slices of the three dimensional MHD turbulence [27]. The cases in which 
the scaling exponents fall in between SL1 and SL2 may be some mixture of the two physical pictures. 
Also, in the global reference frame and in the local perpendicular reference frame the MHD models 
are similar. In these two cases it seems that the dissipation structures are similar for incompressible 
and compressible MHD turbulence. 

HD and MHD turbulence are similar over high order longitudinal scaling exponents measured in 
the global reference frame as well as for motions perpendicular to the local mean field. In the local 
magnetic system the motions parallel to the local mean field for both MHD models are different from 
hydrodynamic motions. 



7 Fast, Slow, and Alfvenic Modes 

The compressible MHD equations support three types of linear waves or modes known as Alfvenic, 
fast, and slow. In the incompressible limit the fast mode does not exist. The velocity fields can be 
decomposed into these three components and scaling exponents for each component can be obtained, 
as explained in [280. Figure [12] gives a pictorial representation of the decomposition method [7j. Figures 
[TBI through [15] contain decomposed scaling exponents and are organized as follows (similar to Figures 

1 Note however that this decomposition method is an approximate one. It is well known that MHD allows 
for discontinuities and shocks. In the vicinity of these structures the wave decomposition method fails. We 
expect that these structures are more pronounced for higher sonic Mach number fluids. Therefore, a proper 
care should be taken for high Mach number fluids (see [7] for an evaluation of the decomposition method 
for a high Mach number case). In our current case, the Mach number is of order unity and, therefore, the 
discontinuous structures should not be very important. 
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[9] through [TTjl . In Figure [13] the global reference frame is used, in Figure [14] the local parallel reference 
frame is used, and in Figure [15] the local perpendicular reference frame is used. In each Figure the 
(a) is the Alfvenic mode, the (b) is the fast mode (containing only the CMHD model), and (c) is the 
slow mode. For the IMHD model the field of the slow mode was calculated by subtracting the Alfvenic 
field from the total velocity field. Each plot in Figures [13] through [15] contains the scaling exponents 
derived from the longitudinal and transverse structure functions in order to expose the differences. For 
the Alfven mode in Figure [T4l given that in the direction parallel to B ext there is no Alfven mode, the 
longitudinal components have no physical meaning and are not included. 

Beginning with Figure [13] in the global reference frame, for the Alfven mode the longitudinal 
components of both models follow SL1 while the transverse follow closer to SL2. For the the fast mode 
the CMHD longitudinal and transverse components are almost indistinguishable and follow SL2. For 
the slow mode the transverse components of both models are more SL2-like. 

In Figure [TT] in the local parallel reference frame, for the Alfven mode the CMHD transverse 
component follows SL2 quite well while the longitudinal component falls between SL1 and SL2. For 
the fast mode the CMHD transverse component falls well below SL2 while the longitudinal component 
falls between SL1 and SL2. For the slow mode the IMHD longitudinal component falls between SL1 
and SL2 while the CMHD longitudinal component follows SL2 quite well. The transverse components 
of both models fall well below SL2 and for the IMHD model the curve actually becomes flat. 

In Figure [15] in the local perpendicular reference frame, for the Alfven mode all cases follow SL1 
except the IMHD transverse component which falls between SL1 and SL2. For the CMHD model the 
longitudinal and transverse components are indistinguishable. For the fast mode both components of 
the CMHD model follow closer to SL2, the transverse component more so. For the slow mode the IMHD 
longitudinal component follows closer to SL1 while the transverse component follows SL2 fairly well. 
For the CMHD model this is the only case where the transverse component is significantly greater than 
the longitudinal component. The transverse component follows SL2 while the longitudinal component 
falls well below SL2. 

There are three cases in which scaling exponents fall well below SL2, namely the CMHD longitu- 
dinal component of C/o C s '!' (p) > the transverse component for both models of Cl°ou>' (p) > an< ^ the CMHD 
longitudinal component of CJio'w (p)> suggesting that intermittency is greater in these cases. 

Also, Cl{p) > Ct(p) for all cases except for the CMHD model of Cllo^ip)- This is not reflected in 
Figure [TT1 where for the CMHD model (l(p) > Ct(p)- This suggests that the amplitude of the slow 
mode for this model is small compared to the amplitude of the Alfvenic and fast modes. The two cases 
where Cl(p) = (t(p) are for the CMHD models of (j l ° st (p) and C^/'^O 5 )- 

In agreement with the above-mentioned property, Cl(p) > Ct(p), it seems that the longitudinal 
components are more likely to be SLl-like while the transverse components are more likely to be 
SL2-like. This is the case for Caifv^P) ano - ^slow which i s m direct agreement with the conspicuous 
cases mentioned in the discussion of Figures [3] through [TT] Combining the results of Figures [5] and [T5] 
it seems that in the global reference frame the Alfvenic mode is mostly responsible for the fact that 
the longitudinal components of both MHD models follow SL1. For the IMHD model both the Alfvenic 
and slow modes seem to contribute significantly to the fact that the transverse components appear 
to follow SL2. Combining the results of Figures [TTJ and Qj}] it seems that in the local perpendicular 
reference frame once again the Alfvenic mode is mostly responsible for the fact that the longitudinal 
components of both MHD models follow SL1. For the IMHD model in this case the slow mode seems 
to contribute most to the fact that the transverse components appear to follow SL2. 

It is also true that for the decomposed scaling exponents of the two MHD models the longitu- 
dinal components are more alike than the transverse components, in every case except for CJi^w (p)- 
This makes sense in the local parallel reference frame where there is a bias in that the longitudinal 
components are chosen in the direction of the local mean field while the transverse components are 
chosen randomly to be perpendicular to the longitudinal components, but this observation also holds 
in the global reference frame where there is no bias because the longitudinal components are chosen 
randomly. If the dissipation structures are similar for compressible and incompressible MHD then it is 
possible that the longitudinal scaling exponents would do a better job of probing this. 
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8 Discussion and Summary 

In this paper turbulence in an incompressible fluid with and without a magnetic field as well as 
moderately compressible MHD turbulence are considered. We summarize our findings as follows: 

1. Substantial similarity was shown between the spectra for all three models. 

2. The PDF's of the velocity components perpendicular to the external magnetic field arc IHD-likc 
while those parallel to the field have a significantly smaller range of velocities. The PDF's of the 
transverse velocity increments for the MHD models decay slower than the IHD model. 

3. The structure functions of the MHD models in the global frame and in the plane perpendicular to 
the local mean field are IHD-like in that the values of the transverse structure functions are greater 
than those of the longitudinal. In these two frames the shapes of the structure functions of the 
IMHD model are more IHD-like than the CMHD model. 

4. The similarity of HD, IMHD, and CMHD turbulence persists over high order longitudinal structure 
function scaling exponents measured in the global reference frame as well as for motions perpen- 
dicular to the local mean field. In these two frames the longitudinal scaling exponents of both 
MHD models seem to follow the predictions for one dimensional IHD-like dissipations structures 
while the transverse scaling exponents of the IMHD model seems to follow the predictions for two 
dimensional IMHD-likc dissipation structures. In the local magnetic system the motions parallel to 
the local mean field for both MHD models are different from hydrodynamic motions. 

5. In the global reference frame and for motions perpendicular to the local mean field the Alfvenic 
mode is mostly responsible for the fact that the longitudinal components of both MHD models follow 
the IHD model and the IHD theoretical dissipation structure predictions. In the global reference 
frame for the IMHD model both the Alfvenic and slow modes seem to contribute significantly 
to the fact that the transverse components follow the IMHD theoretical predictions. For motions 
perpendicular to the magnetic field for the IMHD model the slow mode seems to contribute most 
to the fact that the transverse component seems to follow the IMHD theoretical predictions. It 
is also true that for the decomposed scaling exponents of the two MHD models the longitudinal 
components are more alike than the transverse components, in every case except for the slow mode 
in the plane perpendicular to the local mean field. 
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10° 




Fig. 1 Kinetic energy spectra for the three models as well as Kolmogorov's power-law dependence in the 
inertial range for comparison. 




Fig. 2 PDF's of the velocity field parallel (x) and perpendicular (y & z) to B ex t for the three models. 
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Fig. 3 PDF's, normalized to match their maxima, of the velocity field parallel and perpendicular to B ex t for 
the one time step of the CMHD simulation used throughout the paper and the average over 10 time steps of a 
256 resolution simulation with error bars showing expected random fluctuations. 
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Fig. 4 PDF's of the transverse increments of the x component of the velocity field taken along y and z averaged 
together for the three models. Note that the mean magnetic field is parallel to the x-axis in the MHD cases. 
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Fig. 5 PDF's of the transverse increments of the y component of the velocity field taken along x and z 
presented separately for the three models. 
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Local Mean B Field Local Mean B Field 




Fig. 6 For (a) points are chosen along the local mean field, hereafter referred to as local parallel. For (b) points 
are chosen in the plane perpendicular to the local mean field so that both the longitudinal and the transverse 
structure functions lie in this plane, hereafter referred to as local perpendicular. 
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Fig. 7 Second order structure functions of the velocity field for the three models. The top row are longitudinal, 
the middle row are transverse, and the bottom row are total structure functions. The first column is in the 
global reference frame, the middle row is the local parallel frame, and the bottom row is the local perpendicular 
frame. In the local frames the HD structure functions from the global frame are shown for comparison. 
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Fig. 8 The sixth order longitudinal structure function as a function of the third order longitudinal structure 
function both calculated in the local parallel frame for the three models, al ong with their values of scaling 
exponents obtained from a minimum chi-square linear fit. Compare to Figure [lOl (a). 




Fig. 9 Scaling exponents, normalized to the third order, as a function of the order of the structure functions 
in the global reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, (b) - 
transverse, and (c) - total). 
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Fig. 10 Scaling exponents, normalized to the third order, as a lunction ol the order ol the structure lunctions 
in the local parallel reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, (b) 
- transverse, and (c) - total). For the IHD model the global scaling exponents are shown for comparison. 




Fig. 11 Scaling exponents, normalized to the third order, as a function of the order of the structure functions 
in the local perpendicular reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, 
(b) - transverse, and (c) - total). For the IHD model the global scaling exponents are shown for comparison. 
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Fig. 12 Separation method. We separate Alfven, slow, and fast modes in Fourier space by projecting the 
velocity Fourier component v k onto bases £a, £«, and £/, respectively. Note that £a = —0- Slow basis £ s and 
fast basis lie in the plane defined by Bo and k. Slow basis £ s lies between —8 and ky. Fast basis lies 
between k and kx (from Cho & Lazarian 2003). 




Fig. 13 Decomposed scaling exponents, normalized to the third order, as a function of the order of the 
structure functions of the velocity field in the global reference frame for the IMHD and CMHD models as well 
as K41, SL1, and SL2 (top - Alfvenic, middle - fast, and bottom - slow). 
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Fig. 14 Decomposed scaling exponents, normalized to the third order, as a function of the order of the 
structure functions in the local parallel reference frame for the IMHD and CMHD models as well as K41, SL1, 
and SL2 (top - Alfvcnic, middle - fast, and bottom - slow). 




Fig. 15 Decomposed scaling exponents, normalized to the third order, as a function of the order of the 
structure functions in the local perpendicular reference frame for the IMHD and CMHD models as well as K41, 
SL1, and SL2 (top - Alfvenic, middle - fast, and bottom - slow). 



